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ABSTRACT 

Observations of the optical polarization of NGC 4151 in 1997-2003 show variations of an order of 
magnitude in the polarized flux while the polarization position angle remains constant. The amplitude 
of variability of the polarized flux is comparable to the amplitude of variability of the total U-band 
flux, except that the polarized flux follows the total flux with a lag of 8 ± 3 days. The time lag and the 
constancy of the position angle strongly favor a scattering origin for the variable polarization rather 
than a non-thermal synchrotron origin. The orientation of the position angle of the polarized flux 
(parallel to the radio axis) and the size of the lag imply that the polarization arises from electron 
scattering in a flattened region within the low-ionization component of the broad-line-region. Polar- 
ization from dust scattering in the equatorial torus is ruled out as the source of the lag in polarized 
flux because it would produce a larger lag and polarization perpendicular to the radio axis. We note 
a long-term change in the percentage polarization at similar total flux levels and we attribute this to 
a change in the number of scatterers on a timescale of years. 
Subject headings: galaxies: active — galaxies:quasars:general — polarization 



1. INTRODUCTION 

Shk lovskiH (jl953) pointed out that synchrotron ra- 
diation should be highly polarized, and the detection 
in 1953 of polarization of the Crab Nebula demon- 
strated the synchrotron na ture of its continuum ra- 
diation (jDombrovskiil |1958[ ). It was therefore nat- 
ural to search for polari zation of the optica l emis - 
sion of Seyfert galaxies ( pibai fc Sha khovskovl Il966f l. 
Optical polarization of NGC 41 5 1 was discovered by 
Babadzhanyants & Ha gen-Thornl ([1969) who also re- 
ported that the polarization was va riable. The polar- 
ization was studied in more detail by iKruszewskl (|1971l ) 
who found that the polarization in the V band varied by 
a factor of four over a nine-month period. The detec- 
tion of variable polarization was taken as evidence of a 
synchrotron origin of the continuum. 

The discovery of polarization of emission lines 
(| Angel et all 119761) in NGC 1068 showed that the po- 
larization of type-2 AGNs was due to scattering and 
not synchrotron emission. The origin of the polariza- 
tion in type-1 AGNs, h owever, has remained ambiguous. 
IStockman et al.l (|1979t ) discovered that the optical po- 
larization of radio-loud AGNs was parallel to the radio 
jet axis. They considered this to be either the result of 
intrinsic polarization of optical synchrotron emission, or 
the result of scattering from a flattened distrib ution of 
scatters. Subsequently lAntonuccll (|1982t Il983| ) showed 
that, in general, polarization tends to be parallel to the 
radio axis in type-1 AGNs and parallel in type-2 AGNs. 
Since the polarization of broad lines is roughly similar 
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to th e polarization of the continuum ([Goodrich fc Milled 
1994) this showed that this polarization was due to scat- 
tering. Nonetheless because of the variability of opti- 
cal polarization it has often been assumed, by analogy 
with blazars (where there is no doubt that the optical 
polarization variations have a non-thermal origin) that 
synchrotron radiation also contributes to optical polar- 
ization in non-blazar AGNs (e.g., Giannuzzo & Salvati 
1993). In this paper we attempt to resolve the question 
of the nature of the optical polarization by investigating 
the relationship between the variability of the polarized 
flux and the total flux. 

2. OBSERVATIONS 

NGC 4151 was observed through a 15-arcsecond diam- 
eter aperture with the UBVRI Do uble I mage Chopping 
Photometer - Polarimctcr (Piiro ial fT988h on the 1.25 -m 
AZT-11 telescope of the Crimean Astrophysical Observa- 
tory. The polarimeter gives simultaneous measurement 
of the linear polarization in the standard Johnson bands 
by using dichroic filters. The resulting passbands have 
effective wavelengths of 3600, 4400, 5400, 6900, and 8300 
A which are are close to the standard Johnson UBRVI 
system. Observations of NGC 4151 were made on 58 
nights between May 1997 and May 2003. About 30 
separate measurements were made per night. Errors in 
the flux and polarized flux were estimated by comparing 
measurements taken within ±5 nights of each other. 

3. RESULTS 

We show the variability of t he total ?7-band flux and B- 
band p olarized flux in Fig. 1. iMerkulova fc Shakhovskovl 
( 2006) have already discussed the variability of the to- 
tal flux and polarized flux in the five passbands. As 
is well known, flux variations in the optical passbands 
are nearly simultaneous in AGNs. The [/-band flux 
gives the best indication of the variability of the total 
flux from the AGN since it has minimal host galaxy 
contamination. The polarization might also be ex- 
pected to be strongest towards shorter wavelengths, but 
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Fig. 1. — Variations in total C/-band flux (top curve) and B- 
band polarized flux (lower curve). The polarized flux has been 
multiplied by a factor of 30 for plotting convenience. The two 
curves are fourth-order polynomial fits through the data. 

iMerkulova fe Shakhovskovl (|2006f ) find that the polariza- 
tion variability of NGC 4151 peaks in the B band and 
shows a significant drop in the U band. iKishimoto et al.l 
(|2004f ) find from their spectropolarimetry that the po- 
larized flux spectrum of AGNs commonly peaks in the 
region of the B band and drops off in the U band. They 
attribute this to a Balmer absorption edge in the accre- 
tion disk emission. Because the polarization variability 
of NGC 4151 is at a maximum in the B band and because 
the signal-to-noise ratio is also highest in the B band we 
compare the B-band polari metric variability with the U- 
band total flux variability. ICrenshaw et alJ (| 19961 ) found 
that the UV continuum variations in NGC 4151 are si- 
multaneous to withi n ~ ±0.15 days across the UV and 
lEdelson et alJ (j 19961 ) get an upper limit of ~ 1 d to the 
lag of the A5100 continuum with respect to the con- 
tinuum at A1275. Wav elength-dependent l ags are seen 
at longer wavelengths (Scr geev et alJ 120051 ) . but these 
could be due t o contamination by emission from hot dust 
(Gaskcll 2007). We thus do not expect a significant lag 
of the B band relative to the U band. 

As can be seen from Fig. 1, during the period 
1997-2003 the [/-band flux showed long-term vari- 
ability of a factor of six, and the £?-band polarized 
flux showed long-term va r iabilit y of a factor of ten. 
IMerkulova fe Shakhovskovl (|2006h find that the position 
angle (PA) of the variable polarization remained nearly 
constant at 92 ° ± 1?6 which is consistent w ith earlier 
measurements (|Antonuccil [l983t iMartell I1998D a nd with 
being parallel to the radio axis (|Antonuccilll983l ). 

4. ANALYSIS 

4.1. Long-Term Variability 

During the period 1997 to 2003 NGC 4151 declined 
from the end of a high sta te and passe d through a low 
minimum level of activity (|Lvutvl 120061) . There is some 
evidence for a slow long-term change in the degree of 
polarization over the seven-year period. This is most ob- 
vious if we compare years where the mean t/-band con- 
tinuum level was approximately the same. As can be 
seen from Fig. 1, the polarized flux around MJD 2400 is 
more than a factor of two lower than that around MJD 
1300 when the total fluxes are comparable. Similarly the 
polarized flux is also lower around MJD 2750 compared 
with MJD 1200. This change in the percentage polariza- 
tion after the photometric minimum around MJD 2000 



Fig. 2. — Cross-correlation function for the polarized flux and 
total B-band flux for the entire six-year data set. A positive lag 
corresponds to the polarized flux following the total flux. 

suggests a possible reduction on a dynamical timescale 
in the number of scatterers after the low state. Alter- 
natively, the continuum emission could be intrinsically 
anisotropic and the anisotropy could be varying from 
year to year. Unfortunately we did not have good polari- 
metric coverage in the year of the photometric minimum. 

4.2. Short-Term Variability 

In order to search for polarization reverberation on 
short timescales we have removed the long-term trends 
by subtracting out fourth-order polynomials (shown in 
Fig. 1) from the total flux and polarized flux time series. 
The two time series were cross-correlated usi ng the inter- 
polation method of Gaskcll & S parkd (1 19861) . Deta ils of 
the method are given in iGaskell fc Petersonl (119871) . We 
show the resulting cross-correlation function (CCF) for 
the entire data set in Fig. 2. 

The centroid of the CCF calculated from the first mo- 
ments above 0.7 times the peak correlation gives a lag 
of 8 days. In Table 1 we also give the considerably less 
certain lags from the centroids of the CCFs for the in- 
dividual observing seasons. The median of the lags for 
the individual years is consistent with the lag from the 
whole of the seven o bserving seasons. 

As is well known (|Gaskell fc Petersonl Il987f) the dis- 
tribution of errors in the peaks in cross-correlation func- 
tions is non-Gaussian with a tail extendi ng to high errors. 
This can be seen in the simulations of iMaoz fe Netzen 
(1989), IWhite fc Petersonl (| 19941) . and iPeterson et all 
(1998). We estimated the error in the lag of the po- 
larized flu x in three different ways. First we used the 
formula of IGaskell fc Petersonl (|1987l ). This gave an er- 
ror of ±3.5 d for the lag derived from the whole data 
set. We also performed Monte Carlo simulations where 
we added noise to the observations equal to the obser- 
vational errors. This gave an error in the lag of ±2.5d. 
Finally, where possible we calculated the lags for the in- 
dividual years. They are again the first moment lags 
calculated above 0.7 of the peak in the CCF. These lags 
are shown in Table 1 w i th th e approximate errors from 
the IGaskell fc Petersonl (|1987l ) formula. As can be seen, 
even though the errors are substantially larger because of 
the small number of observations, the median lag for the 
individual years (8.5 d) is close to the lag we obtained 
from all years. The median of the absolute values of the 
differences of the individual lags for each year from the 
lag for all years is 1.5 d. These three separate error es- 
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TABLE 1 

Centroids of Lags for Individual Years 



Year Lag (days) 



1997 7 ±4 

1998 43 

1999 7 ±7 

2000 10 ±5 
2001 

2002 14 ±8 

2003 7 ±5 
All years 8 



timates suggest that the error in the lag for the whole 
sample is of the order of ±3 d. 

5. DISCUSSION 

A ~ 8 d lag of the polarized flux has no natural expla- 
nation if the polarization has a synchrotron origin. How- 
ever, a lag can naturally be explai ned as a light-travel- 
time delay of scattered radiation (|Giannuzzo fc Salvatil 
1993). As noted above, scattering off both electrons and 
dust has long been considered as a cause of polariza- 
tion in AGNs. Electron scattering could take place in 
clouds along the axis of symme try or in a flattened equa- 
torial distribution of electrons (jSmith et al. 2005). Dust 
scattering in the torus will also be responsible for po- 
larization. To investigate which process dominates in 
producing the polarized flux we modeled polarization 
reverberation in NGC 4151 with the Monte Carlo ra- 
diative transfer code STOK ES ((G oosmann fc Gaskelll 
l2007t[Goosmann et al.l l2007). In order to model polariza- 
tion variability we modified STOKES to time-tag each 
photon going through the program. 

5.1. Polarization from the Dusty Torus 

For the period 1990-1998 lOknvanskij et al.l (|1999l ) find 
a delay of 35 ± 8 days for the if -band relative to the U 
band. This includes the peak of activity at the start of 
our polarimetri c monitoring. For th e low state during 
our monitoring iMinezaki et al.l (|2004h find a similar de- 
lay of 48 ± 2.5 days. The K band emission comes from 
the hottest dust, and lags given by the cross correla- 
tion meth od are biased towards material at the small- 
est radii (Gas kell fc Sparkelll986f ) so the A-band lag is 
giving the inner radiu s of the dust torus. From the 
lOknvanskii et al.l {l999) and lMinezaki et alj 1)20041 ) mea- 
surements we can conclude the inner radius of the dusty 
torus during our polarimetric monitoring was ~ 40 light 
days. This strongly rules out dust scattering being the 
c ause of the polariza tion reverberation we detect. 

lOknvanskijl (|1993f ) did determine a K-hand lag dur- 
ing the low state of NGC 4151 in the 1970s of 
18 ± 6 days which is marginally consistent with the 
polarization reve r berati on delay we measure, but as 
Oknvansk ii et al.l (12006D po i nt ou t, their observations 
and those of IMinezaki et al.l (|2004D show that, as would 
be expected, the high state of NGC 4151 at the start of 
our polarimetric monitoring destroyed the dust close to 
the black hole, and the A"-band observations show that 
the dust had still not reappeared at smaller radii within 
several years of the high state. 

We modeled dust scattering off a variety of optically- 
thick torus geometries. We considered cylindrical tori 



and tori with elliptical cross sections. In all cases we set 
the inner radii to be 40 lt-d as indicated by the IR rever- 
beration mapping. Various opening and viewing angles 
were modeled. The resulting first moment delays in the 
polarized flux were ~ 43 days (i.e., slightly greater than 
the inner radius), and varied only by 10-15% with chang- 
ing geometries and viewing angles. These modeled lags 
are clearly inconsistent with the ~ 8 d lag we find from 
the observations. The polarization of the dust model av- 
eraged over the lag was perpendicular to the radio axis 
for all the type-1 viewing angles, which is inconsistent 
with the observed position angle. We thus believe that 
scattering from the torus is not the cause of the lag in the 
polarized flux. For more extensive discussion of the effect 
of the geometry and viewi ng angle on the degree and d i- 
rection of polarization see lGoosmann fc G askell (2007). 

5.2. Polarization from a Flattened Electron Scattering 

Region 

Scattering off a polar distribution of electrons 
(|Giannuzzo fc Salvatilll993l ) also produces the wrong po- 
larization angle, but a flattened electron distribution pro- 
duces polarization parallel to the radio axis (see discus- 
sion in Goosmann & Gaskell 2007). We therefore used 
STOKES to model the polarization lag from cylindrical 
electron-scattering disks. We found that for disks which 
were optically thick to electron scattering the lag was 
equal to the inner radius. Thus if the electron scattering 
region is optically thick in NGC 4151, it has an inner 
radius of ~ 8 lt-d. 

The lag we find is com parable to the size of the BLR. 
iGaskell fc Sparkd (| 19861) obtained a size of 5 ± 2 lt-d 
for the radius of the C IV emitt ing region of th e BLR 
in NGC 4151 during 1978-1980. iMetzroth et "all IpOOfih 
obtained radii of 3.4 ± 1.3d for 1988 and 3.3 ± 0.9d 
for 1991. They a l so ob tained identical radii for He II. 
IGaskell fc Sparkel (|1986l ) estimated the size of the region 
emitting H/3 and H7 to be ~ 6 lt-d durin g 1980 -1981. 
More recent observations by iBentz et al.1 (|2006f ) give 
6.6 ± 1 lt-d for 2005. Note again that the responsivity- 
weighted radii given by reverberation mapping are biased 
towards the inner radii of the emitting region. 

While our observed polarization lag is in good agree- 
ment with the observed radii of the C IV and H/3 emis- 
sion in NGC 4151, the observed inner radius for a given 
line reflects the radial variation of ionization and there is 
almost certainly BLR gas inside that radius. If the den- 
sity of gas in our STOKES models is somewhat lower 
than a canonical BLR density (i.e., 10 8 cm~ 3 rather than 
10 10 cm -3 ), then the model is no longer optically thick, 
the mean free path to electron scattering becomes signif- 
icant, and we find that the lag is greater than the inner 
radius of our hypothetical disk. For example, a disk with 
a density of 10 8 cm~ 3 would give a mean free path of ~ 6 
lt-d. It is therefore possible to reproduce the observed 
polarization lag without having to have an artificial hole 
in the middle of the distribution of electrons. We give a 
more extensive discussion of the results of our modeling 
elsewhere (Shoji, Gaskell, & Goosmann, in preparation). 

Interestingly, the radius we obtain for the electron- 
scattering disk (i.e., towards the outer edge of the BLR) 
is in good agree ment with the locat ion deduced quite in- 
dependently by iSmith et alj (|2005l ) from modeling the 
change in polarization with velocity across the profiles of 
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broad emission lines. 

6. CONCLUSIONS 

The polarized flux of NGC 4151 appears to lag the 
unpolarizcd flux by ~ 8 d which is comparable to the 
light-crossing time of the broad-line region. We inter- 
pret this as the result of the extra light-travel time of 
the scattered photons. Dust in the torus is ruled out 
as the source of polarization by the angle of polariza- 
tion and because the observed lag is too short. Both 
the direction of polarization and the lag can naturally 
be explained instead by electron scattering in a flattened 
region of similar size to the low-ionization BLR. We have 
also detected a possible change in the percentage polar- 



ization of NGC 4151 at similar flux levels on a timescale 
of years, and we attribute this to long-term changes in 
the number of scattering electrons. 
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